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Abstract.  It  has  been  developed  a generalization  of  the  Kohn-Luttinger  envelope-function 
method  that  is  applicable  for  description  of  the  electron  and  hole  states  in  many-layer  (001)  het- 
erostructures, composed  of  related  lattice-matched  III -V  semiconductors,  with  atomically  abrupt 
heterointerfaces.  It  was  shown  that  additional  contributions  to  the  standard  one-band  effective- 
mass  equations  may  be  classified  with  powers  of  the  parameter  kamax  <C  1,  where  l/k  is  the 
characteristic  size  of  the  envelope  function,  and  the  length  alnax  is  of  the  order  of  the  lattice 
constant.  It  was  formulated  a hierarchy  scheme  for  the  effective-mass  equations,  the  «th  level 
of  which  accounts  for  taking  into  consideration  all  corrections  up  to  (few)".  Zero  level  of  the 
hierarchy  corresponds  to  the  standard  effective-mass  equations  with  position-independent  effec- 
tive mass.  On  the  first  level  of  the  hierarchy  each  heterointerface  gives  an  additional  ^-function 
contribution  to  the  potential  energy.  Only  on  the  second  level  the  position-dependent  effective 
mass  appears  as  well  as  corrections  for  the  weak  non-parabolicity  of  the  spectrum  and  spin- 
orbit  interface  interaction.  At  higher  levels  of  the  hierarchy  non-local  contributions  appear,  and  a 
one-band  differential  effective-mass  equation  does  not  exist. 


Currently  the  question  on  applicability  of  the  envelope-function  (EF)  method  of  Lut- 
tinger  and  Kohn  [1]  and  the  effective-mass  (EM)  equations  for  real  nanostructures  is 
discussed  intensively  [2],  There  are  several  problems  associated  with  the  envelope- 
function  method  for  heterostructures.  The  first  one  is  whether  the  method  can  realty 
be  applicable  when  studied  are  the  structures  whose  potential  changes  considerably  on 
the  lattice  constant  scale.  As  the  matter  of  fact  to  researchers’  attention  the  question 
has  been  overtopped  with  the  second  problem  concerning  the  correct  form  of  the  ef- 
fective kinetic  energy  operator  (KEO)  in  the  effective-mass  method,  which  is  just  the 
one-band  EF  approximation,  when  EM  parameters  are  position-dependent.  To  avoid 
this  second  problem  some  authors  prefer  using  complicated  many-band  EF  schemes 
even  when  the  simple  one-band  approximation  would  provide  them  with  at  least  quali- 
tatively correct  solution.  Eliminating  the  KEO  problem  this  way  one  still  has  a trouble 
with  the  first  question.  The  main  goal  of  the  work  is  to  derive  simple  EM  equations 
for  heterostructures  taking  into  account  abruptness  of  the  heterostructure  potential  as 
well  as  position-dependent  EM.  Considered  are  the  (001)  heterostructures  composed 
of  related  (band  offsets  are  small  as  compared  to  the  band  gaps)  lattice-matched  1I1-V 
semiconductors.  The  presented  below  result  is  a kind  of  recipe,  but  the  recipe  based 
on  a rigorous  formalism  of  the  Luttinger-Kohn  EF  method  [2], 

Generally,  one-band  EF  equations  (an  equation  for  simple  band  and  a set  of  equations 
for  degenerate  band)  are  obtained  via  a perturbation  theory  procedure.  So  there  must 
be  small  parameters.  When  one  deals  with  shallow  states,  the  typical  energy  of  which 
E is  small  as  compared  to  the  typical  energy  separation  from  other  bands  Eg,  such  a 

small  parameter  is  Je/E%.  The  small  parameter  may  be  put  down  as  Xk;  here  k is 
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the  characteristic  wave  number  of  the  state  in  question,  and  A = 2TiP/ (moEg),  where 
P is  the  typical  interband  matrix  element  of  momentum  and  mo  is  the  free  electron 
mass.  When  the  proper  band  edge  EM  m is  mainly  formed  with  the  k • p interaction, 

we  may  adopt  A ~ h/  JmE%.  Another  small  parameter  emerges  when  we  face  abrupt 

changes  of  the  heterostructure  potential.  The  parameter  is  ak,  where  a is  the  lattice 
constant.  It  is  important  that  for  a good  deal  of  materials  we  have  a ~ A (e.g.  for 
GaAs  A«6  A).  This  is  the  reason  why  we  should  retain  terms  to  the  same  order  in 
these  two  parameters  in  the  one-band  EF  equations.  Now  we  are  about  to  discuss  the 
obtained  EM  equations.  The  equations  are  presented  the  way  they  govern  electron  and 
hole  states  in  a (001)  symmetric  quantum  well  of  width  L in  the  absence  of  external 
potentials,  and  L = Na/2,  where  A f is  an  integer. 

In  the  usual  (bulk)  EM  equation  [1]  for  electrons 

[T  + V(z)\Fc=EFc, 

the  KEO  and  potential  energy  have  standard  form: 

2 

Tbuik  = Fbulk(z)  = A Ec{0(-z)  + &(z-L)},  (1) 

where  m is  the  bulk  EM,  A Ec  is  the  conduction  band  offset,  Fc  is  the  conduction 
band  EF.  The  function  Fbuik  describes  the  rectangular  shape  of  the  conduction  band 
QW.  Allowing  for  k • p interaction  with  remote  bands  as  well  as  terms  arisen  due  to 
abrupt  change  of  the  potential  at  the  heterointerfaces  give  an  opportunity  to  classify 
these  contributions  with  powers  of  the  parameter  kamax  -C  1,  where  amax  = max(A,  a). 
We  present  the  hierarchy  scheme  for  the  EM  equations,  where  the  nth  level  of  the 
hierarchy  accounts  for  taking  into  consideration  all  corrections  up  to  (^amax)"- 

Conduction  band  QW  states  near  T point 

Oe)  Zero  level  of  the  hierarchy  for  electrons.  The  effective  Hamiltonian  has  “bulk”  form 
(1),  that  is 

f(0)  = 7U;  K(0)(z)  = Fbulk(z). 

We  emphasize  that  at  this  (the  main!)  level  of  the  hierarchy  m(z)  = const. 

le)  First  level  of  the  hierarchy.  As  taken  into  account  are  terms  having  ah  smallness, 
KEO  does  not  change,  and  in  the  potential  energy  there  appear  b-function  corrections: 

T(1)  = 7bu,k;  E(I)(z)  = Ebuik(z)  + d\{5  (z)  + b (z  — /,)}. 

The  appearance  of  the  two  b-functions  at  the  heterointerfaces  is  due  to  abrupt  change 
of  the  potential  and  corrections  arisen  as  we  describe  the  real  (non-abrupt)  form  of  the 
heterostructure  potential  with  the  Heaviside  step-functions.  The  parameter  d\  depends 
both  on  bulk  properties  of  the  structure  components  and  the  microscopic  form  of  the 
transition  region  of  the  heterojunction  (HJ),  and  generally  differs  for  each  different 
heterointerface. 

2e)  Second  level  of  the  hierarchy.  If  we  wish  to  take  into  account  corrections  of 
the  order  of  (ok)2,  (A k)2  and  (ak)(Xk),  we  should  deal  with  the  equation  where  the 
concept  of  the  position-dependent  EM  appears: 

T{2)  = \ma  (z)  pm/j  (z)  p ma  (z)  + a0p4  + (30  (p2p2  + p2p2  + p(p?) ; 
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F(2:  = F(1)  + d 2 [p  x n]  • <j{8  (z)  — 8 (z  — T)} 

Here  m(z)  is  the  position-dependent  band  edge  EM,  a0  and  /3b  are  the  weak  non- 
parabolicity  parameters,  n is  the  unit  Oz-axis  vector,  cr  are  the  Pauli  matrices.  The 
parameters  a (for  (3  we  have  2a  +(3  = — 1 ) and  dj  also  depend  on  bulk  properties  of  the 
materials  of  the  structure  as  well  as  on  the  microscopic  structure  of  the  heterointerface, 
and  should  be  different  for  each  different  HJ  (we  have  a symmetric  structure,  so  the 
parameters  are  the  same  for  both  HJs).  The  parameter  c/2  originates  from  spin-orbit 
interaction,  which  is  taken  into  account  through  the  perturbation  theory,  the  typical 
spin-orbit  interaction  parameter  being  considered  to  be  of  the  order  of  the  band  offsets 
(the  case  takes  place  for  a good  deal  of  heterostructures). 

Now  we  present  the  reasons  why  one  cannot  achieve  better  accuracy  for  the  one- 
band  approximation  in  the  framework  of  differential  equations  of  finite  order.  Originally 
the  many-band  set  of  k • p equations  is  a set  of  k-space  integral  equations,  and  k-space  is 
restricted  with  the  first  Brillouin  zone.  For  shallow  states  one  may  use  an  approximate 
canonical  transformation  of  the  set  (or  another  unitary  perturbation  scheme)  and  move 
onto  the  one-band  EF  approximation.  It  is  important  that  not  all  k belonging  to  the  first 
Brillouin  zone  may  now  correctly  contribute  to  such  an  equation.  The  evidence  comes 
from  the  following.  For  a bulk  material  one  may  use  a series  to  present  the  spectrum 
e(k)  of  the  states  near  T point.  Such  a series  has  a finite  radius  of  convergence  Kq, 
which  is  determined  with  the  strength  of  the  k • p interaction  and  proximity  of  other 
bands.  The  estimation  gives  Kq  ~ 1/A.  This  way  for  the  one-band  EF  equation  k- 
space  is  restricted  with  radius  Kq  rather  than  the  first  Brillouin  zone  boundary.  To 
obtain  differential  equations  in  r-space  one  should  go  beyond  the  k-space  effective 
limits.  The  procedure  would  induce  an  exponentially  small  error  if  either  potential  or 
EF  were  smooth  on  the  lattice  constant  scale.  That  would  mean  the  Fourier  transform 
of  either  potential  or  EF  drops  exponentially  as  k increases.  But  in  our  case  we  have 
abrupt  heterostructure  potential.  One  may  see  from  all  above  that  Fc  is  not  smooth, 
in  particular  its  second  derivative  has  two  discontinuities  with  the  relative  jumps  of 
the  order  of  unit.  This  means  that  if  L is  large  enough,  e.g.  the  conduction  band  EF 
Fourier  transform  Fc( k)  oc  l/(£z)3.  So,  the  principal  error  of  the  one-band  method  is 
(A/ Kq)3  ~ (AA)3.  It  is  the  error  that  does  not  allow  one  to  reach  better  accuracy  than 
that  given  with  the  second  level  equation. 


Hole  QW  states  near  F point 

The  only  difference  between  the  EM  method  for  conduction  and  valence  bands  lies 
in  more  complicated  character  of  the  EF  equations  for  the  latter  case.  This  way  zero 
hierarchy  level  of  accuracy  of  the  one-band  EM  method  is  plain.  As  for  the  second 
one,  the  proper  equation  has  so  many  parameters,  some  of  which  depend  on  specific 
microscopic  details  of  heterointerfaces,  that  in  practice  it  would  be  useless,  and  similar 
equation  for  simple  band  could  have  probably  satisfied  the  academic  interest.  So,  we 
present  the  first  grade.  In  this  approximation  the  Hamiltonian  for  valence  band  states  is 
a sum  of  the  standard  (bulk)  6x6  KEO  matrix,  and  the  potential  energy  6x6  matrix. 
It  is  the  latter  that  is  presented  below.  We  use  the  basis  {| J,jz}}  of  eigenfunctions  of 
the  total  angular  momentum  J and  his  projection  jz,  which  are  linear  combinations  of 
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the  F 15  valence  band  edge  Bloch  functions  and  spin.  We  adopt  the  following  ordering: 
13  3 \ 13  3 \ 13  1\  13  1\  II  1\  |1  1' 


1 2’  2/  ’ 1 2’  2/’  1 2 ’ 2 / ’ 1 2’  2/’  1 2’  2/  ’ 1 2’  2, 

and  the  phase  agreement  is  the  same  as  in  [1],  So,  the  potential  energy  matrix  is 


V = 


Frsl 

F0V„ 

y'v/21/0t(jv 


V0CTy  -iV2V0CTy 

Frgl  0 


(2) 


0 


VTl\ 


Here 

Fr8  = EVi  + A(/r8{0  (-z)  + 0 (z  - L)}  + Xl  (z)  + <5  (z  - T)}, 

Fr,  = ETl  + Af/r,{©  (— z)  + © (z  - L)}  + X2{«5  (z)  + <5  (z  - L)}, 

Vo  = X3{<5  (z)  — 8 (z  — L)}. 

The  notations  are  the  following:  Ey,  and  Ey.  are  the  Tg  and  IT  band  edge  energies 
of  the  well  material,  Af/r8  and  At/r,  are  the  proper  band  offsets.  Two  parameters  Xi 
and  X2  have  the  same  origin  as  d\  for  the  conduction  band  case,  and  the  parameter  xi 
appears  only  as  we  correctly  process  the  abrupt  changes  of  the  heterostructure  potential 
(that  means  the  parameter  vanishes  for  smooth  potentials).  'Hie  last  parameter  defines 
in  particular  the  strength  of  the  zone  center  mixing  of  light  and  heavy  holes  [3]. 

In  conclusion  we  presented  three  grades  of  accuracy  that  can  be  achieved  in  the 
one-band  EF  approximation,  and  this  way  formulated  three  hierarchy  levels  for  the 
EM  equations.  There  are  effects  that  cannot  be  described  using  the  zero  level  of  the 
hierarchy,  e.g.  zone  center  light  hole-heavy  hole  mixing  in  the  valence  band.  There 
are  also  effects  that  can  be  dealt  with  only  using  the  second-level  equations  of  the 
hierarchy,  e.g.  spin  splitting  of  the  subbands  when  longitudinal  electric  field  is  applied 
(when  the  presence  of  the  terms  like  the  one  proportional  to  c/2  is  crucial),  or  some 
optical  transitions. 

The  work  was  supported  by  RFBR  (project  96-02-18811),  RFBR-INTAS  (95-0849), 
the  Federal  Programs  “Physics  of  Solid  State  Nanostructures”  (96-1019),  and  “Surface 
Atomic  Structures”  (95-3.1). 
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